Universal contact of strongly interacting fermions at finite temperatures 
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The recently discovered universal thermodynamic behaviour of dilute, strongly interacting Fermi 
gases also implies a universal structure in the many-body pair-correlation function at short distances, 
as quantified by the contact X. Here we theoretically calculate the temperature dependence of this 
universal contact for a Fermi gas in free space and in a harmonic trap. At high temperatures above 
the Fermi degeneracy temperature, T > Tf, we obtain a reliable non-perturbative quantum virial 
expansion up to third order. At low temperatures we compare different approximate strong coupling 
theories. These make different predictions, which need to be tested either by future experiments or 
advanced quantum Monte Carlo simulations. We conjecture that in the universal unitarity limit, the 
contact or correlation decreases monotonically with increasing temperature, unless the temperature 
is significantly lower than the critical temperature, T<T C ~ 0.2Tf- We also discuss briefly how to 
measure the universal contact either in homogeneous or harmonically trapped Fermi gases. 

PACS numbers: 03.75.Hh, 03.75.Ss, 05.30.Fk 



I. INTRODUCTION 

Understanding strongly interacting fermions is one of 
the most challenging problems in present-day physics []| . 
While the behavior of interacting fermions can be under- 
stood in some well-defined regions of parameter space, 
it remains elusive in the correlated strongly interacting 
regime A better understanding of strongly inter- 

acting fermions has wide-ranging implications for sys- 
tems such as quark matter in neutron stars and high- 
temperature superconductors. An important generic idea 
in this field is fermionic universality 0, H[ : all strongly 
interacting, dilute Fermi gases should behave identically, 
depending only on a scaling factor equal to the average 
particle separation, but not on the details of the interac- 
tion. 

The recent realization of broad collisional (Feshbach) 
resonances with an interaction potential range r$ small 
relative to the inter-particle spacing in ultracold atomic 
Fermi gases provides a highly controlled environment 
for studying the general problem of strongly interact- 
ing fermions [H, [J. By applying an external magnetic 
field across the resonance, the interparticle interaction 
can be accurately tuned from weak to infinitely strong 
[6]. This has led to the observation of crossover from 
Bardeen-Cooper-Schrieffer (BCS) superfluids to Bose- 
Einstein condensations (BEC) [mI- 

The most strongly interacting regime lies at the reso- 
nance, where the s-wave scattering length a s diverges and 
the two-body scattering amplitude reaches its maximum 
value allowed by quantum mechanics, i.e., it becomes 
unitarity limited [4] . A number of properties in the uni- 
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tarity limit have been characterized. In particular the 
universal thermodynamic behavior [H, [Iol - [l5| has been 
experimentally observed, and agrees with approximate 
analytic theories and Monte Carlo simulations. However, 
it is still nontrivial to calculate these universal properties 
quantitatively, as there is no small interaction parameter 
in a perturbation theory expansion. 

In 2008, Tan gave new insight into this difficult prob- 
lem by deriving a set of relations which link the short- 
distance, large momentum correlations to the bulk ther- 
modynamic properties of a fermio n sy stem with short- 
range interparticle interactions |16l - l3ll |. These relations 
generically apply to any dilute Fermi gas with an inter- 
particle spacing much larger than Tq. All of these Tan re- 
lations are exact in the limit of a vanishingly small range 
of the interaction potential, so that kpro — > 0, where 
hp is the Fermi momentum. They are connected by a 
single coefficient X, referred to as the integrated contact 
intensity or "contact". For instance, the pair correlation 
function is predicted to diverge as Xj (I67rr 2 ) at short dis- 
tance r(3> rn) — > [l6| and the momentum distribution 
will fall off as I/A; 4 at large momentum fc(<C 1/Vo) [l7| . 
The original definition of the contact given by Tan is, 



X = lim 

k— >oc 



P<r(k), 



(1) 



where p a (k) is the momentum distribution in one spin 
component. The contact X is an extensive quantity and 
has the unit of Nkp, where N is the total number of 
atoms of the system. For a homogeneous Fermi gas, it 
is also convenient to use a contact intensity, C = X/V, 
where V is the volume of the system. 

The immediate significance of the Tan relations can 
be seen most clearly from Tan's adiabatic sweep theorem 
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which states that the rate of adiabatic change of the to- 
tal energy with respect to the inverse scattering length is 
proportional to the contact. While it is surprising that 
the short-range, high-energy correlations can be related 
so seamlessly to the low-energy equation of state, the 
Tan relations are closely related to the Feynman-Hellman 
theorem[32j]. This method has also been used to calculate 
correlations in one-dimensional Bose gases at finite tem- 
perature from thermodynamic properties [33], with exper- 
imental verification 34 1. 

In the fermionic strongly interacting regime, universal 
thermodynamics implies a universal contact. The deter- 
mination of the contact therefore provides a very impor- 
tant means of characterizing the many-body properties 
and phases of strongly interacting fermions, complemen- 
tary to existing measurements of the equation of state. 
As it is related to the derivative of the observed energy, 
the contact is a more sensitive test of theoretical pre- 
dictions than direct thermodynamic measurements. We 
note that some of the Tan relations have been confirmed 
experimentally at JILA in the USA [35| and at Swinburne 
University of Technology in Australia (3(| . 
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Figure 1. (Color online) The experimental measured contact 
for a trapped interacting Fermi gas (symbols) is shown as a 
function of the dimensionless coupling constant l/(fc_pa s ). It 
is compared with a zero-temperature perturbation prediction 
(solid line) [36| and a finite-temperature second-order virial 
expansion calculation at 0.5Tf (dashed line) [2Ej |. The mea- 
surements include Bragg spectroscopy (empty squares) [3^ |. 
momentum distributions (solid or empty circles) |35l |. rf spec- 
troscopy (stars) [351 ] . and molecular fraction (triangles) [38l |. 

Experimentally, the contact is obtained by a number 
of methods, including measuring: 

• the molecular fraction [37l . f38j | , 

• the momentum distribution tail (35| . 

• the radio-frequency (rf ) spectroscopy signal at large 
frequencies [35l . [39j , or 

• the spin-antiparallel static structure factor at large 



Figure 1 summarizes all the measured contact data to 
date for a trapped interacting Fermi gas (symbols). The 
experimental results are compared with our earlier pre- 
dictions of a zero-temperature result from an approxi- 
mate Gaussian pair fluctuation theory (solid line) [36j 
and from a finite-temperature calculation using the quan- 
tum virial expansion (dashed line) [25| . The experimen- 
tal data of course are not taken exactly at zero temper- 
ature. Accurate determination of temperature is diffi- 
cult due to the absence of an efficient thermometry. A 
typical experimental estimate gives T/Tp ~ 0.10 — 0.20. 
The data in Fig. 1 agree reasonably well with the zero- 
temperature curve, although near the unitarity limit they 
seem to lie systematically lower, possibly due to the 
nonzero temperature in these experiments. A similar 
theoretical calculation at zero temperature has been re- 
ported by Werner and collaborators [38| by interpolating 
the known perturbation results in the BEC and BCS lim- 
its. 

In this paper, following the pioneering studies by Yu, 
Palestini and co-workers, we systematically investigate 
the temperature dependence of the contact for the most 
interesting case of a unitary Fermi gas. Yu et al. es- 
timated the finite-temperature behavior of the universal 
contact (40| . based on an assumption of phonon excita- 
tions at low temperatures and a second-order quantum 
virial expansion at high temperatures; while Palestini et 
al. presented ab-initio results using a non-self-consistent 
T-matrix approximation (4l| . Both studies predicted 
a pronounced maximum in the contact at finite tem- 
peratures of about 0.6T^ and 0.2Tp, respectively. The 
present investigation, however, calls into question this 
prediction of a temperature-dependent peak in the uni- 
versal contact. Our main results are summarized below, 
following a two-fold motivation. 

Firstly, we give a high-temperature quantum virial ex- 
pansion study of the contact up to third order. This 
provides us a general idea of how the virial expansion 
converges with decreasing temperature, and therefore en- 
ables us to estimate the temperature window for its ap- 
plicability. We find that the virial expansion of the con- 
tact is quantitatively accurate down to the Fermi degen- 
erate temperature Tp and 0.5Tp for a homogeneous and 
trapped Fermi gas at unitarity, respectively. We show 
that, in exact analogy with the virial expansion of ther- 
modynamic potential [42| . the contact at unitarity can 
also be expanded in terms of universal, temperature in- 
dependent coefficients, c„ = dAb n /d(X/a s ), referred to 
as the contact virial coefficients. Here, Ab n is the n-th 
virial coefficient and A is the thermal wavelength. This 
is a direct consequence of fermionic universality. We pre- 
dict that for a uniform unitary Fermi gas, the second 
and third contact coefficients are given by, respectively, 
C2.00 = 1/tt — 0.318 and 03.00 ~ —0.141. The subscript 
"00" stands for the unitarity limit with a s = ±00. 

Secondly, we are interested in the low-temperature be- 
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havior of the universal contact. In the strongly inter- 
acting regime, as there is no controllable small inter- 
action parameter, a comparative study using different 
strong coupling theories is necessary |43l445j . We find 
that the theoretical descriptions of the finite-temperature 
contact from different strong coupling theories show con- 
siderable discrepancies near the critical temperature T c 
(~ 0.2Tp) for the onset of superfluidity. In particular, 
the enhancement (peak structure) of the contact near 
T c , found earlier in a non-self-consistent T-matrix calcu- 
lation [ill] , is completely absent in other strong coupling 
theories. Therefore, we conjecture that the universal con- 
tact decreases monotonically with increasing temperature, 
unless the temperature is much smaller than the critical 
temperature. This discrepancy remains to be resolved by 
future experiments and advanced quantum Monte Carlo 
simulations [46l-{48j|, exploring the critical temperature 
regime close to resonance. 

Our paper is organized as follows: in Sec. II we present 
the quantum virial expansion for the contact and calcu- 
late the universal contact coefficient at unitarity, for both 
trapped and homogeneous Fermi gases. We then discuss 
how to determine the contact from strong-coupling theo- 
ries in Sec. Ill, using different means via the tail of mo- 
mentum distribution, Tan's adiabatic sweep relation and 
pressure relation. Theoretical predictions for the uniform 
universal contact from different strong-coupling theories 
are compared with each other and with accurate virial 
expansion results in the high-temperature regime. In the 
next section (Sec. IV), we describe how to calculate the 
trapped universal contact from the homogeneous data 
and report the results of different strong-coupling theo- 
ries and virial expansion. In Sec. V, we discuss briefly 
the measurement of the temperature dependence of the 
universal contact in a homogeneous Fermi gas. A sum- 
mary and outlook is given in Sec. VI. 



all the quantum states of n-particles with a proper sym- 
metry. The virial coefficient b„, can be calculated from 
the cluster partition functions [42l. [52j. for instance, 



and 



b 2 = Q2/Q1 - Qi/2 



b 3 = Q3/Q1 -Q2 + Qi/3, etc. 



(4) 



(5) 



In practice, it is often convenient to separate out a 
background, non-interacting thermodynamic potential, 



-ksTQ^z + b { 2 1] z 2 



b { n ] z n 4 



Here, 



we have used the superscript "I" to indicate the non- 
interacting systems. The thermodynamic potential of in- 
teracting Fermi gases can then be rewritten as, 



fi = - k B TQ 1 [Ab 2 z 2 + ■■■ + Ab n z n 



(6) 



where Ab n = b n — 6„ '. For a homogeneous Fermi gas, 
Qi = 2V/X 3 with volume V and thermal wavelength A = 
fiTrtf/imkBT)] 1 / 2 . 

The essential assumption of the quantum virial expan- 
sion is that the above expansion of n - converges 
for temperatures down to T c , the critical temperature 
for the onset of superfluidity. We note that the fugac- 
ity in the expansion may be larger than unity at these 
low temperatures. However, convergence is still possible 
given small enough virial coefficients Ab n . For a homo- 
geneous unitary Fermi gas, the second virial coefficient 
Ab2.oo — 1/V2 was known 70 years ago (55|. The third 
virial coefficient, A&3 j00 ~ —0.35501298, was recently cal- 
culated by the present authors [42] and confirmed experi- 
mentally in an accurate thermodynamic measurement by 
Nascimbene and co-workers (l3| . 

The virial expansion of the contact follows directly 
from an alternative representation of Tan's adiabatic 
sweep theorem in the grand- canonical ensemble, 



II. QUANTUM VIRIAL EXPANSION OF THE 
CONTACT 



The quantum virial expansion has proved to be an ef- 
ficient method for studying the high-temperature prop- 
erties of ultracold atomic Fermi gases |42l . |49to4| . This 
method utilizes the fact that in the high temperature 
limit, the chemical potential \x diverges to —00 and the 
fugacity z = exp(/i/fceT) < 1 is a well-defined small ex- 
pansion parameter. Thus, one may expand any physical 
quantities of interest in powers of fugacity, no matter how 
strong the interaction is. 

For the grand thermodynamic potential of a two- 
component Fermi gas, the expansion is given by (42l l52l|. 



= -k B TQi [z + b 2 z 2 



H h b n z n + > 



(3) 



where Q n — Tr„[exp(— "H/fcsT)] is the partition function 
of a cluster containing n particles and b n is the n-th virial 
expansion coefficient. The trace Tr n takes into account 
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(7) 



This is simply because the adiabatic sweep theorem im- 
plies the first law of thermodynamics, 

AE = n 2 I/(4 7 rm)A(-a7 1 ) + TAS + fiAN , (8) 

which can alternatively be written as 

Afi = ft 2 I/(47rm)A(-a7 1 )-S'AT-iVA/i. (9) 

Therefore, using Eq. ([6]) we immediately obtain a quan- 
tum virial expansion for the contact: 
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+ c n z n + 



(10) 



where we have defined the dimensionless contact coeffi- 
cient, c„ = dAb n /d(X/a s ). For a homogeneous system, 
we sometime use the contact intensity, C = I/V. 
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In general, the contact coefficient should be a func- 
tion of X/a s and hence is temperature dependent. In the 
unitarity limit where \/a s = 0, however, we anticipate 
a constant, universal contact coefficient, similar to the 
universal virial coefficient A&„ j00 |42l. |49| . This is a man- 
ifestation of fermionic universality, shared by all systems 
of strongly interacting fermions 0, 



A. Universal relation between the homogeneous 
and trapped contact coefficients 

In exact analogy with the virial coefficient [42l |. 
fermionic universality leads to a very simple relation be- 
tween the trapped and homogeneous contact coefficients 
at unitarity. Let us consider the contact of a harmon- 
ically trapped Fermi gas with the trapping potential, 
V T (r) = muj 2 r 2 /2 or V T (r) = muj 2 ± (x 2 +y 2 ) /2+mu 2 z 2 /2 
with uj = (cj^cj^) 1 / 3 . In the thermodynamic limit of 
lu — > 0, we may use the local density approximation 
and neglect the discrete energy levels. The whole Fermi 
system is treated as many cells with a local chemi- 
cal potential /x(r) = /i — Vt{y) and a local fugacity 
z(r) = e ^ r )/ fcsT = zexp(-V T /k B T). Due to the con- 
stant contact coefficients, the spatial integration in the 
total contact It = J dr[C(r)] can be easily performed. 
Here, we use the subscript "T" to indicate explicitly the 
trapped system. We find that, 



Airm 



kBTQi iT X [c 2 ,oo,tz 2 + c 3i00i t-z 3 H ] 



(11) 

where the trapped contact coefficient is given by a uni- 
versal relation, 



Cn,oo,T 



7,3/ 2 ' 



(12) 



and Qi,t = 2(fcsT) 3 / (hui) 3 is the single-particle parti- 
tion function in harmonic traps and in the local density 
approximation. 

In the following, using the known solution of two- and 
three-fermion problems, we calculate the universal second 
and third contact coefficients, in both homogeneous and 
trapped configurations. 



B. Second universal contact coefficient 

The second contact coefficient of a homogeneous inter- 
acting Fermi gas can be obtained from the well-known 
phase- shift expression for the second virial coefficient 

SHI, 



^h. = e -Ei/(k B T) 
V2 V 



i 



d6 Q 
~dk 



(13) 



Here, is the energy of the i-th bound state of the 
two-body attractive interaction and So(k) is the s-wave 



scattering phase shift, given by k cot <5o — — l/a s +rok 2 /2. 
To obtain the derivative dAb2/d(\/a s ), let us choose the 
BCS side with a s < 0, on which the two-body bound 
state is absent. As the range of the interaction potential 
ro is vanishingly small so that fc^T <C f?j(mr^) and 
k <C l/?'o, we find that ddo(k)/dk ~ — a s /[l + k 2 ^] and, 



A6 2 (a s < 0) 



V2 



dt 



1 



1 + t 2 



-A 2 t 2 /(2™') 



(14) 



Near to the resonance, it is readily seen that A62 — 
1/V2 + A/(?Ta s ), giving rise to a homogeneous contact 
coefficient, 



1 



(>2, 



(15) 



To calculate the trapped second contact coefficient, we 
consider the second virial sufficient in an isotropic har- 
monic trap, which is given by (42I IH^. 



1 



A&9.T = - 



2 ^ 



el,n/k B T 



} „/k B T 



(16) 



Here. 



„ = (2v n + 3/2)Huj is the n-th relative energy of 
two fermions with unlike spins and v n satisfies the secu- 
lar equation, 2T{—v n )/T(—v n — 1/2) = d/a Sl with F the 
gamma function and d = \j2h/ (mui) the characteristic 
length scale of the harmonic trap. The non-interacting 
relative energy is e^] n = (2n + 3/2)hui (n = 0, 1, 2, ...) 
and in the unitarity limit, the solution of v n is known 
analytically, v n .oo = n — 1/2. It is easy to show that, 



de r 



d(X/a s 



\/a B =0 



HujdT(n+ 1/2) 
?tA n! 



Thus, we find that in the unitarity limit, 



= d - ^ r O + l/ 2 ) r -(2n+im 
2ir\ 



C2,oo,T-— 



(17) 



(18) 



n=0 



where u) = huj/(kBT) <C 1 is the reduced trapping fre- 
quency. The sum over n can be exactly performed, lead- 
ing to, 

-1 1/2 



1 



C2,oo,T — 



2^2-k 

1 

2V2n 



2uj 



(19) 



The leading term in the above equation is universal, sat- 
isfying the universal relation Eq. (|12p . The second term 
(oc uj 2 ) is non- universal and is caused by the length scale 
of the harmonic trap [42j . It represents the finite-size cor- 
rection to the local density approximation that we have 
adopted above. This correction however is extremely 
small if u> — > 0, as anticipated. At the Fermi degen- 
erate temperature, Tp = (3_/V) 1//3 fiu;/fcB, we find that 
Co 2 cx N~ 2 / z ~ 10~ 4 , for the typical number of atoms 
N r~j 10 5 in the current experiments |2j. 
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C. Third universal contact coefficient 

The determination of the third contact coefficient is 
more cumbersome. In the homogeneous case, it is ex- 
tremely difficult to use the conventional method em- 
ploying the three-particle scattering matrix to obtain 
the third virial coefficient. Therefore, it is more prac- 
tical to determine firstly the trapped contact coefficient 
C3.00.T, and then to use the universal relations at low 
trap frequency to obtain the homogeneous result, 03.00 = 
3v3c3 )00 ,T- We note that this method relies on exact 
three-body solutions to the energy eigenvalues in a trap. 
These are known analytically for all eigenstates, and can 
be summed numerically to any required accuracy. 

An estimate of C3 j00j t can already be obtained by the 
known results of Ab 3>F as a function of the coupling con- 
stant l/k F a s at different temperatures T/T F and at a 
small S> ~ 0.1 (see, for example, Fig. 36 in Ref. (4^|). 
simply because, 

= 1 g>A6 3 ,t f^T~ dAb 3 , T 
C3 ' T ~ k F Xd(l/k F a s ) V 4nT F d(l/k F a s )' [ ' 

We find that the coefficient c^ t00t T at resonance is indeed 
nearly temperature independent and estimate from the 
slope of Abs^r that, C3 iOC) r (estimate) ~ —0.0265 at w ~ 
0.1. 

An accurate determination of c^ t00t T requires a sys- 
tematic extrapolation to the limit of Hi = 0. For this pur- 
pose, we calculate numerically the derivative C3 >00) t(<&) = 
[dAb3 y T/d(X/a s )]x/a s =o as a function of ui. Using the 
small Co data, a numerical extrapolation to u) = gives 
rise to the trapped third virial contact coefficient, 

c 3 ,oo,t =s -0.0271. (21) 

Thus, we obtain immediately from the universal relation, 
Eq. (fT2")l . the homogeneous third contact coefficient, 

c 3 ,oo = -0.141. (22) 

D. Large-T contact: the homogeneous case 

We are now ready to calculate the universal contact in 
the high temperature regime. For a homogeneous Fermi 
system, the single-particle partition function Q\ = 2V/A 3 
and the dimensionless contact T/ (Nk F ) is given by, 

(23) 

Here N = pV is the total number of atoms with the 
homogeneous density p. 

The fugacity z is determined by the number equation 

m, 

p = ~p {1) (z) + [2A6 2)00 z 2 + 3A6 3i0O z 3 + •••], (24) 



where we have defined a dimensionless density p = 
pA 3 /2 = [4/(3V^)](T F /T) 3/2 and the density of a non- 
interacting Fermi gas as 

/■OO 

p (1) {z)= (2/V5F) / dtVt/(l + z- 1 ^) . (25) 

In practice, for a given fugacity, we calculate the di- 
mensionless density using Eq. ([24)) and hence the re- 
duced temperature T/T F . The dimensionless contact is 
then obtained from Eq. (|2"3")h as a function of T/T F or 
the inverse fugacity z _1 . 
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Figure 2. (Color online) Universal contact of a homogeneous 
unitary Fermi gas at high temperatures, as predicted by the 
quantum virial expansion method up to the second order 
(dashed line) and the third order (solid line). Dashed verti- 
cal line indicates the Fermi degenerate temperature T F . The 
inset shows the contact as a function of the inverse fugacity. 

Figure 2 reports the temperature (main figure) or fu- 
gacity (inset) dependence of the homogeneous contact in 
the unitarity limit, calculated by virial expanding to the 
second order (dashed line) or third order (solid line) . The 
close agreement between the second and third predictions 
strongly indicates that the virial expansion works quan- 
titatively well down to the Fermi degenerate temperature 
T F , as indicated by the vertical dashed line. At sufficient 
high temperatures, where 

z^p=[A/(i^)]{T F /Tfl\ (26) 

the leading temperature dependence of the contact is 
given by, 

-4- (T » Tp) = — ( — ) , (27) 
Nk F K ' 3 \T F ) ' v ' 

as predicted by Yu and co-workers [4(| ■ We note however 
that the prefactor there is smaller by a factor of 47r 2 , due 
to a different definition for the contact. 
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E. Large-T contact: the trapped case 

For a trapped Fermi gas at unitarity, the dimensionless 
contact can be written as, 



It „ 3/2 / T \ 7 ^ 2 r 2 q 

— = 24. [-1 [c 2iOC , T z + c 3)00 , T z 



(28) 



The number equation takes the form [52j, 
~p T = ~p ( t ] 0) + [2A6 2 . oc , T z 2 + 3A63 !0O , T z 3 + ■■■], (29) 

where p T = (N /2)(hwf / \k B Tf = (T F /T) 3 /6 and the 
density of a non-interacting trapped Fermi gas p^ (z) = 
(1/2) Jg 00 dtt 2 /(l + z- 1 e t ). The trapped virial coefficient 
is given by A&„ j00j t = A&„ !00 /n 3 / 2 . In analogy with the 
homogeneous case, for a given fugacity we determine the 
reduced temperature T/Tp from the number equation 
and then calculate the trapped contact using Eq. 
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Figure 3. (Color online) Universal contact of a trapped uni- 
tary Fermi gas at high temperatures, obtained by expanding 
the virial series to the second order (dashed line) and the third 
order (solid line). The inset shows the contact as a function 
of the inverse fugacity. 

Figure 3 presents the virial expansion prediction for the 
trapped universal contact, expanding up to the second or- 
der (dashed line) or third order (solid line). Amazingly, 
because of the factor of n~ 3 / 2 reduction for the n-th con- 
tact coefficient in harmonic traps, the convergence of the 
expansion is much improved. The expansion now seems 
to be quantitatively reliable down to 0.5T F . The asymp- 
totic behavior of the contact at very high temperatures 
can be determined by setting 



z — Pt — (T F /T) 3 /6 . 



We find that, 



T 



-5/2 



(30) 



(31) 



in agreement with our previous result from the large mo- 
mentum behavior of static structure factor [25]. Thus, 
the contact in harmonic traps decays at high tempera- 
tures much faster than in homogeneous space, due to the 
reduction of the peak density at the trap center at high 
temperatures. 



III. LOW-T HOMOGENEOUS CONTACT FROM 
STRONG COUPLING THEORIES 

At low temperatures, we have to resort to strong cou- 
pling theories to determine the contact. In the absence 
of a controllable small interaction parameter, however, 
there is no known a priori theoretical justification for 
which strong coupling theory is the most appropriate 
[45l | . Here we consider three different theories within 
the many-body T-matrix approximation: the Nozieres 
and Schmitt-Rink (NSR) Gaussian pair fluctuation the- 
ory [7l. [56l - f59l | . non-self-consistent GnGo theory [o0l-l62l|. 
and self-consistent GG theory [63-65]. 



A. Brief review of T-matrix approximations 

The T-matrix approximation is also referred to as lad- 
der approximation, which involves infinite set of Feyn- 
man diagrams — the ladder sum in the particle-particle 
channel. It is now generally accepted that this ladder sum 
is necessary for accounting for strong pair fluctuations in 
the strongly interacting re gim e and is the leading class 
of all sets of diagrams (66l l67l| . The need for including 
infinite sets of Feynman diagrams is self-evident, since we 
are dealing with a strongly interacting theory in which 
there is no small coupling constant that would allow a 
truncation of perturbation theory to finite order. 

There are several T-matrix theories, differing in the di- 
agrammatic structure of the T-matrix t{Q), the particle- 
particle propagator x(Q)> an d the self-energy 'S(K). 
The ones we consider here are respectively given by the 
following form [66l , |67| , 



U 



m [i + u x (Q)y 

* = E K g « (*) G - w - K ) 

Z(K) = y* t(Q)G a (Q - K). 



(32) 

(33) 
(34) 



Here and throughout, Q — (q, iu n ), K — (k, iui m ), while 
U~ x = m/ (Anh 2 a s ) — ^ k l/(2ek) is the bare contact in- 
teraction renormalized in terms of the s-wave scattering 
length. We follow the conventional notations and use 
e k = h 2 k 2 /(2m) and ^2 K = k B T^2 m ^2 k , where q and 
k are wave vectors, while v n and oj m are bosonic and 
fermionic Matsubara frequencies. The Greek subscript a 
in the expressions for the propagator \ (Q) an d the self- 
energy £ (K) may either be set to "0", indicating a non- 
interacting Green's function Gq(K), or be absent, indi- 
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eating a full interacting Green's function G (K). The lat- 
ter can be calculated once the self-energy is determined, 
using the Dyson equation, 



G{K) = 



Go (K) 



[1 - Go (K) E (K)] 



(35) 



By taking different values of a, two different T-matrix 
approximations are obtained: either self-consistent (no 
a), or non-selfconsistent (a = 0) . We note that, in the 
superfluid phase the Green's function has to be a 2 x 2 
matrix, accounting for the U(l) symmetry breaking. 

The simplest choice in the T-matrix approximation 
is to take a non-interacting Green function Gq(K) ev- 
erywhere in x (Q) an d E(if). This idea was pioneered 
by Nozieres and Schmitt-Rink (NSR) for an interacting 
Fermi gas in its normal state 0] • In addition, an approx- 
imate Dyson equation for the Green's function, i.e., 



G (K) = G (K) + G (K) E (K) G (K) . 



(36) 



was used. Interestingly, this non-self-consistent NSR ap- 
proach combined with a truncated Dyson equation can be 
reformulated using the functional path-integral language 
for the grand thermodynamic potential [56] , expanded to 
second order to include Gaussian fluctuations. From this 
perspective, the NSR theory is the first step in an order- 
by-order expansion in path-integral fluctuations, which 
has a systematic path to including higher-order terms. 

The NSR theory was extended recently to the broken- 
symmetry superfluid phase by several authors in different 
manners [57I . I59L l69l - l71j , some of which involve assump- 
tions in order to ease the computational workload. A 
full extension following the original NSR approach was 
reported by the present authors in 2006 [57], with the 
use of a mean-field (2x2 matrix) BCS Green's function 
as "Go". Hereafter we shall refer to this extension as a 
Gaussian pair fluctuation theory or the GPF approach. 

The truncation of the Dyson equation, Eq. (f3"fj|) . can 
be avoided, as shown by Strinati [60l 161] and Combescot 
[62| and their co-workers. In particular, in the superfluid 
state the coupled T-matrix equations have been solved 
[6l| . In the following, we shall refer to this non-self- 
consistent T-matrix approximation with no truncation in 
the Dyson equation as the GoGo theory. However, since 
this approach is no longer part of a systematic expansion 
of the path integral, it is an open question as to whether 
it is more or less accurate than the NSR theory. 

More sophisticated strong-coupling theories can also be 
obtained by using the full (dressed) Green function G (K) 
in x (Q) an d E (K). This modification, referred to as the 
self-consistent GG theory, was investigated in detail at 
the BEC-BCS crossover by Haussmann and collabora- 
tors |63l . |65| . above and below the superfluid transition 
temperature. One advantage of the GG approximation 
is that the theory satisfies the so-called $— derivability 
and is thus conserving, although the conservation is at 
the single-particle level only. The GG theory has also 
been intensively studied in the condensed matter physics, 



particularly for high-temperature materials, under the 
name of the fluctuation exchange (FLEX) approxima- 
tion [66|, [67J • Below threshold the GG approximation is 
difficult to treat at unitarity, and generally involves ad- 
ditional assumptions in order to obtain a scale-invariant 
theory. 

It is clear that in both the NSR or Go Go approxi- 
mations one omits infinite classes of diagrams that are 
responsible for multiparticle interactions. The fully self- 
consistent GG theory attempts to correct for this, by 
modifying the single-particle Green function in the lad- 
der diagram. However, the more crucial interaction ver- 
tices remain unchanged. At this stage, there are no gen- 
eral grounds for deciding which strong-coupling theory 
is the most appropriate. Therefore, any theoretical pre- 
dictions produced by different strong coupling theories 
should be treated conservatively and be examined criti- 
cally using more accurate ab-initio quantum Monte Carlo 
simulations or future experiments. We note that, for the 
equation of state of strongly interacting fermions, a com- 
prehensive comparison between different strong-coupling 
theories and most recent ex per imental results was re- 
ported by the present authors [451 ] . In this earlier compar- 
ison we also included intermediate (GGo) theories which 
are partially self-consistent (68|. These are omitted here 
for simplicity. 

In the following, we calculate the universal contact at 
unitarity by using the GPF or NSR theory and compare 
our results with these reported by Strinati and co-workers 
(Go Go) [4l| and with the predictions provided by Hauss- 
mann and co-workers (GG) [65l. |72|. |73| . 



Determination of the contact 




Figure 4. (Color online) Large wave-vector power-law behav- 
ior in the momentum distribution p a (k), obtained using differ- 
ent strong-coupling theories. The temperature T — 0.3Tp is 
larger than the critical temperature T c ~ 0.2Tf- The dashed 
lines indicate the asymptotic value of 3iv 2 k 4 p a {k)/k%, which 
gives rise to the dimensionless contact X/(Nkp). The inset 
shows the plot in the linear scale. 



s 



There are a number of ways to calculate the contact 
from T-matrix theories according to the different Tan 
relations, as follows: 



f 



Following Tan's original definition [16j, the most 
direct way seems to be the use of the large wave- 
vector asymptotic behavior of the fermionic mo- 
mentum distribution, Eq. ([T]). The theoretical con- 
tact from the GoGo theory was obtained in this 
manner by Strinati and co-workers [4l|. In the 
JILA experiment, the contact was extracted by av- 
eraging k A p a {k) over an interval at kc,expt ~ 2fcp 
[35( 1 . Figure 4 shows k 4 p a (k) for a unitary Fermi 
gas in its normal state at T = 0.3Tp, obtained 
by sol ving (iteratively) the coupled T-matrix equa- 
tions [64], Eqs. (52 1) . (P |) . ([5^ 1 . and ([55 ]) . The 
power-law fc" 4 tail is clearly identified since k A p a (k) 
approaches a constant above a characteristic wave- 
vector kc ~ 5k p , which is however larger than 
kc,expt- The small wiggles in Fig. 4 at about lOkp 
are due to the numerical inaccuracy. We note that 
the power-law tail of /c~ 4 in the momentum dis- 
tribution was first pointed out by Haussmann |63| , 
using the self-consistent GG theory for a normal 
interacting Fermi gas. In the GG theory, the tail in 
p a {k) implies that the contact can be directly cal- 
culated from the order parameter A and the vertex 
function IYx, t). 



2 _ 3n 2 [A 2 -r u (x = 0,r = 0-)] 



Nk F 



(37) 



A brief explanation of this expression is given in 
the Appendix A. We also refer to Sec. IID of Ref. 
[72| for more details. 

2. Another way to calculate the contact is to use Tan's 
pressure relation [13], 



PV 



-E 



h 2 T 
127rma s ' 



(38) 



which results directly from the adiabatic sweep 
relation. As discussed in Ref. [Tt]], this is ev- 
ident if we write the total energy in the form, 
E = Ne F f E [l/(k F a s ),S/(Nk B )}, where f E is a 
dimensionless function. The adiabatic sweep re- 
lation leads to, h 2 I/(ATrm) = — NeFf' E /kF , where 
the derivative is taken with respect to the variable 
l/(kpa s ). We then find that, 



A direct application of the pressure relation in the 
unitarity limit is prohibited as a^ 1 = and there- 
fore the dependence on the contact is lost. How- 
ever, we can use the pressure relation to calculate 
the contact around the unitarity and then do a 
smooth interpolation to the limit of l/(kpa s ) = 0. 

3. Finally, the contact can also be calculated from the 
entropy by using, 



dl 
df 



47T77J 



dS 



di-as 1 ) 



(42) 



7> 



which can be obtained directly by taking tempera- 
ture derivative in Eq. |7[). This important relation 
was first shown by Yu and co-workers (40| . Thus, 
the temperature dependence of the contact is deter- 
mined by the variation in the entropy with respect 
to the coupling constant. 



Variation in the entropy 

O Tan's pressure relation 
□ n (k) tail 




1.5 



0.0 



0.5 



T/T 



1.0 



1.5 



Figure 5. (Color online) Temperature dependence of the con- 
tact in the unitarity limit, obtained by using the variation in 
the entropy Eq. (|42[) . Tan's pressure relation Eq. (|38[). and 
the momentum distribution tail Eq. Jl]). Here we use the 
Gaussian pair fluctuation theory [53]. 

Figure 5 reports the universal contact calculated using 
the above mentioned methods, within the GPF (NSR) 
theory. There is an apparently discontinuous behaviour 
at the critical temperature T c , due to the breakdown 
of the NSR approach just above T c [13, liij . We find 
that all three theoretical methods yield nearly the same 
prediction, owing to the consistency in thermodynamic 
relations. The self-consistent GG theory has the same 
advantage of thermodynamic consistency. 



PV 



-V 



-N 



dE 



dV 
127rma s 



f E - Nepf E 



y d_k^ 
dV 



1 



(39) 

(40) 
(41) 



C. Universal homogeneous contact from 
strong-coupling theories 

Figure 6 presents the comparison of different theoret- 
ical predictions for the temperature dependence of the 
contact of a homogeneous Fermi gas at unitarity. The 
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Figure 6. (Color online) Universal contact of a homogeneous 
Fermi gas in the unitarity limit. The predictions from differ- 
ent strong coupling theories are compared with each other. At 
high temperatures, they are compared with the virial expan- 
sion results as well. The symbols at low temperature indicate 
the predictions calculated by using the measured equation of 
state (star) [T3 |. the quantum Monte Carlo result for pair cor- 
relation function or ground state energy (circle) [74l. Tr^| and 
the lattice quantum Monte Carlo simulation at the critical 
temperature (square). In the inset, the temperature depen- 
dence of the contact is shown for a BCS Fermi gas at the 
interaction strength l/{kpa s ) = — 1. The circle at T = 
shows the contact obtained from the quantum Monte Carlo 
result for pair correlation function [75j |. 



X/(Nk F 



Methods 



3.51 ± 0.19 ENS equation of state data (Ref. [14]) 
3.396(9) QMC ground state energy (Ref. [77]) 
3.40 QMC pair correlation function (Ref. [741 1 
3.34 GPF (present work) 
3.02 GG (Ref. [72]) 
3.23 Go Go (Ref. [41]) 

Table I. Zero temperature contact in the unitarity limit. The 
contact T/{Nk F ) ~ 67r£/5 can be calculated from the ground 
state energy of an interacting Fermi gas near unitarity [l7| : 

E/(Ne F ) = (3/5) [£ - (/{k F a s ) H ], where £ and £ are two 

universal parameters. The contact may also be determined 
from the tail of spin-antiparallel static structure factor [25l ]. 
Sti-il) — Z/(ANq). Here, we compare different results ob- 
tained from the three strong-coupling theories with the pre- 
dictions by using £, either measured at ENS (£ = 0.93(5)) 
P3] or extracted from quantum Monte simulation for ground 
state energy (£ = 0.901(2)) [77j or for the static structure 
factor (or pair correlation function (£ = 0.90)) [74l ]. 



zero temperature results are listed in the Table I, as well 
as the contact calculated by using the low temperature 
experimental data of the equation of state [l4| and by 
using the quantum Monte Carlo data for pair correla- 
tion function [74| or ground state energy [77|. The ex- 



perimental data of course is not at zero temperature, 
as this is not experimentally achievable. In these mea- 
surements, low temperature thermometry was difficult. 
Similar experiments carried out elsewhere have reported 
lowest achievable temperatures of around T/Tp — 0.05. 

At low temperatures, we observe two distinct predic- 
tions for the behavior of the contact with increasing tem- 
perature: while the GqGq theory predicts an increase of 
the contact and hence a maximum around the critical 
temperature T C) both the GPF and GG theory suggest 
that the contact decreases monotonically as the temper- 
ature increases, unless T -C T c . This qualitative discrep- 
ancy deserves a careful analysis. 

There are known arguments for the presence of a max- 
imum in the contact at low temperature, as follows: 

• The enhancement of the contact just above T c in 
the Go Go theory was interpreted to be related to 
the strengthening of local pairing correlations in the 
absence of long-range order (4l| . 

• As shown by Yu and co-workers j4jj, phonon exci- 
tations at low temperatures will enhance the con- 
tact as T 4 , causing a growth in the contact with 
temperature. 

• The assumption of Fermi-liquid behavior for a 
weakly coupled Fermi gas [40] would also lead to 
a maximum in the contact at T ~ Tp. 

These arguments, however, are not convincing enough to 
conclude there is a real maximum in the low-temperature 
contact in the unitarity limit, as we discuss below. 

The difficulty is that all the strong coupling theories 
are less accurate near criticality, suffering from strong 
pair fluctuations. Thus, the enhancement of the contact 
just above T c in the GqGq theory might be related to 
the breakdown of the approximations in this approach. 
Indeed, we observe a very similar enhancement of the 
contact in the NSR theory. However, we know that the 
NSR approach simply breaks down in the temperature 
region just above criticality. 

On the other hand, while phonon excitations are im- 
portant at low temperatures, our realistic calculation us- 
ing the GPF theory at low temperatures suggests that 
the power-law of T 4 due to phonons is exhibited at 
T < 0.07Tf < T c only, as clearly shown in Fig. la. At 
such low temperatures, the contribution from phonons is 
of relative order 10 -3 . This is so nearly negligible that 
it might be easily be compensated by other excitations 
in this strongly interacting supcrfiuid at unitarity. 

We have also investigated the Fermi-liquid theory of 
the contact for a weak-coupling Fermi gas, as reported in 
Fig. 7 b, where we compare the Fermi-liquid result with 
the NSR prediction for a weakly interacting Fermi gas 
with kpa = —0.5. The essential idea of the Fermi- liquid 
theory is that at low temperatures, the entropy density is 
given by s = m*kpk' 1 B T j (3h 2 ), where the effective mass 
to* = TO,[l+8(71n2-l)/(157r 2 )A:|a2]. Therefore, by using 
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At high temperatures, we find that the NSR and GqGq 
theories predict essentially the same result, as the dif- 
ference between the two theories becomes rather small 
in the high-temperature regime. However, both predic- 
tions lie systematically below the quantum virial expan- 
sion result, indicating the inaccuracies of these theories. 
In contrast, at high temperatures the self-consistent GG 
result approaches the second-order virial expansion pre- 
diction. None of the strong-coupling theories can pro- 
duce correctly the third-order expansion result, since the 
many-body T-matrix approach fails to account for the 
three-particle scattering process. 



Figure 7. (Color online) (a) In the unitarity limit, the low 
temperature contact calculated from the GPF theory (solid 
line) is compared with the T 4 power-law of phonon excitations 
(dot-dashed line). The inset shows the comparison in a larger 
temperature window, (b) The contact for a weakly interacting 
Fermi gas at the coupling constant l/(fci?a a ) = —2. The NSR 
result (solid line) is contrasted with a prediction from the 
Fermi- liquid theory (dashed line). 



Eq. , the derivative of the contact with temperature 
is given by, 



d [1/ (Nk F )} 
d (T/T F ) 

This leads to, 
1 2 



16(71n2 - 1)tt „ .3 T „ 
V 15 (k F a s f — > 0. (43) 



Nk f 



Nk f 



8(71n2-l)7r , t ,/T 
15 {kFas) ll> 



(44) 



where the zero temperature homogeneous contact Iq can 
be calculated using the Lee and Yang's ground state en- 
ergy (up to (kpa s ) 2 ) and has the form 



Nk~ F = 3 {kFas) 



12(11 -21n2), 
1 + — itz L k F a 6 

357T 



(45) 



It is clear from Fig. 7b that the weak- coupling Fermi- 
liquid prediction works at very low temperatures only 
and strongly over-estimates the contact once T > 0.2T F , 
although the theory describes qualitatively the correct 
behavior of an increasing contact up to T ~ T F . Near 
the unitarity limit, the use of Eq. (|44[) and hence its 
support for an enhancement of the contact is certainly 
doubtful. 

In brief, we feel that in the unitarity limit it is more 
reasonable to expect a monotonically decreasing contact 
as the temperature increases. It is interesting to note 
that, for a weaker coupling constant, l/(k F a s ) = —1, 
all the strong-coupling theories give the same qualitative 
behavior of the contact, as shown in the inset of Fig. 6. 
In particular, the contact is predicted to decrease with 
temperature in the superfluid phase. It is natural to as- 
sume that the same decrease would happen as well for a 
unitary Fermi superfluid. 



IV. LOW-T TRAPPED CONTACT FROM 
STRONG COUPLING THEORIES 

Let us now turn to the contact of a unitary Fermi gas in 
a harmonic trap. Under the local density approximation 
this is given by 



drC(r) = (3^ 2 ) 1/3 / dr[C/(/>*F)](r)p 4 / 3 (i 



(46) 

where p(r) is the density distribution and k F (r) = 
[37T 2 p(r)]-V3. 

In the unitarity limit, the local contact [C/(pk F )](r) 
is a function of the density only, through the reduced 
temperature T/T F (r). The density may be expressed 
using p = 8/(3 v 4 : )A" 3 (T F /T) 3 / 2 , so that the trapped 
contact is given by, 



T T = — 

3A 4 



47rr dr 



pk F } 



rp2 



(47) 



Here, for simplicity, we assume an isotropic trap. In the 
local density approximation, it is convenient to use a local 
inverse fugacity C (r) = e~^/ kBT = Co exp[U T (r) /k B T] 
with Co = e~ A1 / fcsT . The integration over the radius r 
can then be converted to an int egration about th e in- 
verse fugacity, by using r = (2fcsT/mo; 2 )£(C) with 
£(C) = ln(C/Co)- I n a harmonic trap, the total num- 
ber of atoms is related to the Fermi temperature by 
N = (l/3)(fcsTp) 3 /(Su;) 3 and the Fermi wave-vector is 
given by k F = y // 2mkBT F /h 2 . Thus, we find that, 



( x 



\Nk t 



32 T 7 / 2 



T 



7/2 



£(c) .(o4 (c) 



Co 



f 2 {c y 



(48) 



where the inputs, the local dimensionless contact I = 
C/{pk F ) and the local reduced temperature T = T/Tp, 
are functions of the inverse fugacity C, and can be de- 
termined using the homogeneous equation of state. The 
temperature T/Tp in the above equation (|48[) is yet to 
be related to the inverse fugacity Co- F° r this purpose, 
we rewrite the number equation N = J drp(r) into the 
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form. 



^ ] = - / d( 



T 



Co 



d( f(C) 3/2 ' 



(49) 



For a given Co, we use Eqs. gHD and (@g) to calculate 
the reduced temperature T/T F and the trapped contact 
[I/(Nk F )]t, respectively, and in turn express the contact 
as a function of the reduced temperature. The detailed 
procedure of numerical calculations can be found in the 
Appendix B. 



(T/(Nk F )) T 


Methods 


3.00 ±0.12 


Swinburne Bragg spectroscopy (Ref. [36J) 


3.40 ±0.18 


ENS equation of state data (Ref. [14]) 


3.26 


GPF (present work) 


3.03 


GG (Ref. [72]) 


3.05 


G Go (Ref. [41]) 



Table II. Zero temperature trapped contact in the unitarity 
limit. We compare different results obtained from the three 
strong-coupling theories with the experimental measurements 
at Swinburne \M and ENS fijl. 



A. Trapped contact at zero temperature 

At zero temperature, the trapped contact at unitarity 
can be found analytically. In this case, the density profile 
is known exactly, 



P (r) = Po 1 



R 2 

T F 



3/2 



(50) 



where the peak density hq and the Thomas-Fermi Rtf 
are respectively given by, 



- £-3/4 



P0 = £ 



(24iV) 
3^ 



1/2 



h 

mu! 



-3/2 



Rtf = £, 



1/4 



(247V) 1 / 6 \ /A 
V muj 



(51) 



(52) 



where £ ~ 0.41(1) is the universal parameter of a unitary 
Fermi gas Because the temperature is zero, the local 
dimensionless contact [C / (N pk)]h om is spatially indepen- 
dent, as it can only depend on a scale factor which has 
already been included in the definition of this quantity. 
Here, for clarity we explicitly indicate the homogeneous 
contact by the suffix "horn". 

Thus, the trapped zero-temperature contact is given 
by, 



1 T 



VMWhom J 



105 \pkF /hom 



Po n TF- 



(53) 
(54) 



Using kp = \j2mksTp /h 2 in harmonic traps with 
k B T F = (SiV) 1 / 3 ^, we find then, 



( 1 



\Nk F 



256 1/4 f_C_\ = 512C 



(55) 



where the two universal parameters £ and £ are related to 
the ground state energy of an interacting Fermi gas near 
unitarity 0, E/(Ne F ) = (3/5)[£ - C/(k F a s ) + ■ ■ ■]. 

Table II lists the zero-temperature contact of a uni- 
tary Fermi in harmonic traps. The theoretical predic- 
tion from different strong-coupling theories are compared 



with the experimental measurements at the Swinburne 
Unive rsity of Technology, Melbourne [36| and at ENS, 
Paris [lj|. We note that, the Swinburne data point is 
determined from spin-antiparallel static structure factor 
measured by Bragg spectroscopy. The lowest tempera- 
ture in the Bragg measurement is about 0.10 ± 0.02Tp, 
which could explain in part the slightly lower value ob- 
tained for the contact. For the ENS data point, we 
have used £ = 0.41(1) and £ = 0.93(5), as determined 
from the experimental data for ground state energy. The 
temperature of this measurement was below the limits 
of the thermometry used, but we can estimate this as 
0.05 ± 0.05Tf from comparisons with related experimen- 
tal measurements. 



Universal trapped contact from strong-coupling 
theories 
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Figure 8. (Color online) Universal contact of a trapped Fermi 
gas in the unitarity limit, calculated from different strong- 
coupling theories as well as from the quantum virial expansion 
method. The symbols indicate the measurements at ENS, 
Paris (star) [3| and at the Swinburne University of Technol- 
ogy, Melbourne (square) (3§ |. 



Figure 8 reports the temperature dependence of the 
contact of a trapped Fermi gas at unitarity, obtained from 
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different strong-coupling theories. We compare the the- 
oretical results with two experimental data points mea- 
sured at low temperatures, while at high temperatures, 
the results are contrasted with the accurate quantum 
virial expansion prediction. 

At low temperatures, T < 0.3Tp, there is a sizable dis- 
crepancy between the predictions from different strong- 
coupling theories. Nevertheless, all the predictions are 
now in qualitative agreement, as the maximum in the 
homogeneous contact found previously in the GqGq the- 
ory is washed out by the trap averaging [4l|. At high 
temperatures, the difference between different theories 
become much smaller. This is in line with our obser- 
vation for the n-th contact coefficient, which receives a 
factor of 7i 3 / 2 reduction in harmonic traps. 

The trapped contact of a unitary Fermi gas can be 
measured readily in experiment, by extending the exist- 
ing low-temperature measurements of either the momen- 
tum distribution |35| or Bragg spectroscopy of a unitary 
Fermi gas [36] to the high-temperature regime. The pre- 
cise determination of the temperature in the unitarity 
limit is a challenging problem [l^ . However, this can 
be overcome by using the known equation of state in the 
unitarity limit [l3l[43ll45| . We anticipate that the quanti- 
tative discrepancy between different strong-coupling the- 
ories shown in Fig. 8 will be clarified by future experi- 
ments. 



V. EXPERIMENTAL MEASUREMENT OF THE 
HOMOGENEOUS CONTACT AT UNITARITY 

It would be very useful to determine the contact of a 
homogeneous unitary Fermi gas, without the complica- 
tion caused by the trap averaging. The result would pro- 
vide a unique opportunity to test quantum many-body 
theories of strongly interacting Fermi gases that we have 
discussed above. However, an in-situ measurement is 
generally required, which is difficult to setup experimen- 
tally. The determination of the universal homogeneous 
contact could be achieved as follows, by: 

• Measurement of homogeneous equation of state 
near unitarity. The Salomon group at ENS, Paris 
has already demonstrated the accurate measure- 
ments of uniform equation of state either in the 
unitarity limit [l3T | or at zero temperature [bj ] . The 
measurement could be extended straightforwardly 
to the general case with arbitrary coupling con- 
stant and temperature, say, for example, the en- 
tropy S(T/Tp,l/kFa)- The homogeneous contact 
can then be determined directly using Eq. P2")h 

• Measurement of the tail of the spatially-resolved 
RF spectrum at unitarity. Tomographic RF spec- 
troscopy was recently developed at MIT to measure 
the local pairing gap of a unitary Fermi gas of 6 Li 
atoms at finite temperatures [7g. Ideally, the lo- 
cal universal contact could be retrieved from the 



tail of RF-spectrum that obeys a w~ 3 / 2 power-law. 
However, the analysis of the tail structure is com- 
plicated for 6 Li atoms due to a strong final state 
effect. It would be desirable to develop a spatially- 
resolved RF-spectroscopy for 40 K atoms, for which 
the final state effect may be negligible. 

Measurement of the spatially-resolved Bragg scat- 
tering spectrum at unitarity Thus, one 
obtains the local spin-antiparallel static structure 
factor at large momentum q, which at unitar- 
ity satisfies the asymptotic behavior of S^i(q) = 
[Z/(7Vfcj?)]/(4(7/fcj?). The local universal contact is 
then determined. 



^ 2 
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Figure 9. (Color online) The universal contact of a homoge- 
neous unitary Fermi gas as a function of the inverse fugacity, 
calculated using the Gaussian pair fluctuation theory, self- 
consistent GG theory and quantum virial expansion. 

All of these measurements are based on the local den- 
sity approximation. The local contact is then to be deter- 
mined as a function of the local chemical potential fj,(r) 
or local inverse fugacity £(r) = exp(— / u(r)//csT). In Fig. 
9, we show the inverse fugacity dependence of the uni- 
versal contact of a unitary Fermi gas, to be confronted 
with future experiments. 



VI. CONCLUSIONS 

To summarize, we have presented a comparative the- 
oretical study of the universal contact for a strongly in- 
teracting Fermi gas at unitarity. The contact at high 
temperature has been accurately determined by using a 
quantum virial expansion method, while at low temper- 
atures, we have employed different strong-coupling the- 
ories to try to estimate the temperature dependence of 
the contact. 

The temperature dependence of the contact near the 
critical temperature, however, remains unresolved, since 
one of the strong-coupling theories, the GqGq theory, 
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predicts an enhancement or a maximum near criticality, 
which is not observed in the other two strong-coupling 
theories. We believe that this enhancement could be due 
to the inaccuracy of the GqGq approach and conjecture 
that the universal contact should decrease monotonically 
with increasing temperature, except possibly at very low 
temperatures T <C T c . 

As there is no solid justification for any of these strong- 
coupling theories, our conjecture should be examined 
critically by future experiments or advanced quantum 
Monte Carlo simulations. Therefore, we have proposed 
several experimental ways to determine the homogeneous 
contact of a unitary Fermi gas, using measurements of 
equation of state, tomographic RF-spectroscopy, and 
spatially-resolved Bragg spectroscopy. All of these mea- 
surements are within reach of current experimental tech- 
niques. 

We note finally that fermionic universality and the 
detailed behavior of a strongly interacting Fermi gas 
near the normal-superfluid transition is of great interest. 
Tan's contact provides an entirely new means to char- 
acterize the universal properties and phases of strongly 
interacting fermions, in addition to the equation of state. 
Our comprehensive study of the universal contact, built 
on the most recent theoretical methods, should provide 
useful insights for future research. 

Note added. Since the submission of our 

manuscript, the finite-temperature contact of a trapped 
unitary Fermi gas was measured at Swinburne [80|. The 
experimental result agrees fairly well the theoretical pre- 
dictions in Fig. 8. A quantum Monte Carlo simulation of 
the finite-temperature contact of a homogeneous unitary 
Fermi gas was also carried out by Drut, Lahde and Ten 

m. 
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Appendix A: Contact from the vertex function 

In the many-body theory, the contact is given by the 
vertex function in the limits of short distance and time. 

For a normal Fermi gas, this can be understood by 
calculating the large-frequency tail of the momentum 
distribution p a (k). For this purpose, we consider the 
diagram contribution to AG(i,t = _ ) = G(x,t = 
Q — ) — GW(i,t = - ), as shown in Fig. 10, where 



-k,-i(0 



k, i co, 




k, ico K 



(AG) 



Figure 10. Diagram contribution to AG(x,t) — G(x,t) — 
G^(x,t) at short distance and time. The shadow region 
indicates the vertex function at x = and r = CP. 



G^(x,t) is the finite-temperature Green function of an 
ideal Fermi gas. In the limit of large momentum k — > oo 
or short distance x — > 0, we may assume that the vertex 
function T(x,t — _ ) is a smooth function and there- 
fore set T(x,t = 0") ~ r(0,0"). According to Fig. 10, 
the finite temperature Green function in the momentum 
space is then given by (p = 0), 



AG(M^m) ~-r(o,o-)— 



i 



i 



(Al) 



where uj m — (2m + VjitkgT is the fermionic Matasubara 
frequency and = H 2 k 2 / (2m). The momentum distri- 
bution Ap a (k) = kBTJ2 m AG {k,iuj m ) takes the form, 



A Pa {k) ~ -r(o,cr)- 



1 



1 



e k /(k B T) + l\ 4e 2> 



(A2) 



(A3) 



Thus, the contact is given by, 1 — —m 2 T(0,0 )/fi 4 or 
1 _ 3tt 2 T(x = 0,r = 0") 



Nkr. 



(A4) 



In the superfluid phase, we should explicitly introduce 
the pairing gap and use a 2 by 2 matrix for the Green 
function and vertex function. This leads to Eq. (|37l) . 



Appendix B: Calculation of the integrals in the local 
density approximation 



The two integrand functions in Eqs. (|48|) and (|49|) . 
which involve the local reduced temperature, are ill- 
behaved at low temperatures. Therefore, we regularize 
the reduced temperature by its non-interacting value and 
define, 



*(C) 



Tig(C)' 



(Bl) 
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where Tig (C) i s the reduced temperature of an ideal, 
non-interacting Fermi gas, defined by: 



3>c(C) = [/(0] 



-2/3 



where: 



(B2) 



respectively. Here, {Ci} (i = 1, iV) are the set of points 
where i (C) has values. We shall let Co run over all the 
points of {d}. To calculate A (Ci), because of the smooth- 
ness of t (C) we approximate it by, 



The coupled equations then take the form, 



(B3) 



i T _ 32 r dy/m, (n ,4/3 fn * (0 , R/n 



Co 



T|_ 32 f,,dy/HO 

T3 ~ 7T ./ ^ dC 

Co 



nC)/(C)[i(C)]" 3/2 , (B5) 



where £ (C) = In(C/Co), and /(C)/* 2 (C) and [t(C)]^ 3/2 
are now smooth functions of C- 

To do the integration, for instance, for the integral in 
the number equation (|B5|) . 



(B6) 



Co 

we use the following discretized version, 

A = E Hd)+ioo, 



(B7) 



where the two types of integrals I (Ci) and are defined 
by, 

Ci+i 



.!(,;)= / dC^2nC)/(C)[t(C)]" 3/2 , (B8) 



Ci 



Aoo = j d ^^^nof(0[t(or 3/2 , (B9) 

c« 



Ci+i 



A (Ci) « [* (C)P 3/2 / d t^r^ e (0 / (0 , (BIO) 

Ci 



[*(0] 



-3/2 



! / (Coe^) , (Bll) 



where [t (CF 3/2 = [t^ 2 (Ci) + i" 3/2 (Ci+i)]A and in 
the second line we have changed to a new variable x = 
\J I (C). For the integral Aoo, where Cjv is sufficiently 
large, we may use the virial expansion result, i.e., 



[*vo (C)]- 3/3 = 1 + A6 3 ,ooC- 3 ^ (m2) 



with Afe 2 ,oo = 1/V2 and A6 3j00 = -0.35510298. Thus, 
we have, 



da;/ Coe x 



-3/2 



(B13) 



The integral in the trapped contact equation (|B4j) can be 
done in a similar way. We note that at high temperatures, 



no 

t 2 (C) 



37T 2 [c 2 ,ooC 2 + C 3 ,aoC~ 



[/(C)] 



4/3 



-, (B14) 



with C2,oo = 1/7T and C3 iQO = —0.141. 
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